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Abstract 

A subtractionless method for solving Fermi surface sheets (FSS), from 
measured n-axis-projected momentum distribution histograms by two- 
dimensional angular correlation of the positron-electron annihilation 
radiation (2D-ACAR) technique, is discussed. The window least squares 
statistical noise smoothing filter described in Adam et aL, NIM A, 337 
(1993) 188, is first refined such that the window free radial parameters 
(WRP) are optimally adapted to the data. In an ideal single crystal, the 
specific jumps induced in the WRP distribution by the existing Fermi 
surface jumps yield straightforward information on the resolved FSS. 
In a real crystal, the smearing of the derived WRP optimal values, which 
originates from positron annihilations with electrons at crystal imper- 
fections, is ruled out by median smoothing of the obtained distribu- 
tion, over symmetry defined stars of bins. The analysis of a gigacount 
2D-ACAR spectrum, measured on the archetypal high-T c compound 
Y Ba2CusOj^s at room temperature, illustrates the method. Both 
electronic FSS, the ridge along TX direction and the pillbox centered 
at the S point of the first Brillouin zone, are resolved. 
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1 Introduction 

While the technique of two-dimensional angular correlation of the positron- 
electron annihilation radiation (2D-ACAR) [p], ||, |J, was successfully used to 
resolve the Fermi surface topology in various classes of materials (see, e.g., f|] 
for a recent review), the attempts to use it in high-T c superconductors met 
huge difficulties stemming from various sources: the occurrence of a fraction 
of Fermi electrons which is significantly lower than in normal metals, the 
layered crystalline structures of the samples, the high proportion of crystal 
imperfections in the single crystals under investigation. 

A straightforward consequence of the occurrence of a layered structure, 
which prevents the uniform distribution of the thermalized positrons inside 
the sample, is the limitation of the usefulness of the 2D-ACAR technique to 
the resolution of the electronic Fermi sheets only. The occurrence of a small 
fraction of Fermi electrons results in very weak Fermi surface jumps in the 
measured momentum density, which can thus be easily missed unless the 
crystal quality is very high and the accumulated statistics is large enough to 
get the small Fermi surface jumps resolved. 

The resolution of the Fermi surface in high-T c superconductors is a chal- 
lenge irrespective of the used technique, however. Each of the techniques 
used so far (2D-ACAR, de Haas-van Alphen, or angle-resolved photoemission 
spectroscopy (see, e.g., J4| and references quoted therein), has been able to 
resolve only part of the existing Fermi surface sheets, such that the final rep- 
resentation of the Fermi surface came from the superposition of essentially 
disjoint pieces of information |J. 

In this paper, we consider the identification of the Fermi surface jumps 
from a 2D-ACAR histogram which records the projection, along a principal 
crystallographic axis n, of the momentum density coming from zero-spin 
positron-electron pairs. Such a histogram will be called thereafter an n-axis- 
projected histogram. 

If the single crystal used in a 2D-ACAR experiment would be perfect, then 
the statistical noise smoothing by means of an appropriate method (see, 
e.g., |0, 0, U), possibly combined with a method of subtracting the radially 
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isotropic component, would resolve the characteristic Fermi surface jumps. 

In fact, there is an important discrepancy between the experimentally 
measured and theoretically computed momentum distributions. In the range 
of low momenta, the latter, which are calculated under the assumption of 
perfect crystalline periodicity, are found to be sensibly smaller than the for- 
mer H . The origins of this discrepancy have been traced back to the positron 
annihilation with electrons at crystal imperfections |5|. A quantitative 
estimate |j] suggested that more than one third of all annihilations belong 
to this category. 

To rule out the spurious effects coming from crystal imperfections, either 
their contribution is to be explicitly subtracted from the data (as done, e.g., 
by the Argonne group JTU], §J), or the smoothing method is to be insensitive 
to the occurrence of crystal imperfections. 

Here we report a method of the second kind. There are two basic ingredi- 
ents of this method: (i) the derivation of optimally adapted to the data radial 
parameters of the constant-weight window least squares (CW-WLS) smoothing 
method reported in || and (ii) removal of crystal imperfection dependence 
of the obtained values by median smoothing over stars of symmetry relating 
the bins of the histogram. 

The occurrence of Fermi surface jumps in the data results in jumps in 
the optimal radial parameters of the CW-WLS filter. The presence of crystal 
imperfections blurs the occurring jumps. The median smoothing restores the 
jumps, if any. 

Essential for the success of such an analysis is the operation of histogram 
redehnition from the laboratory frame (LF), Op^Py p® , (which is defined by 
the setup and within which the data acquisition is performed) to the crystal 
frame (CF) (which is identified with Tp x p y p z , the canonical reference frame of 
the first Brillouin zone of the crystal JTT]| and within which the various steps of 
the off-line analysis of the spectrum are legitimate). Here, this preliminary 
problem is assumed to be solved along the lines described in |]12| Within 
the procedure described in these papers, the assumption (following from the 
experiment) of a momentum projection along a principal crystallographic 
axis is duly checked, such that the expected symmetrization operations are 
justified and hence the inclusion of stars of bins into the smoothing procedure 
of the median method is a valid operation. 

The paper is organized as follows. In Section 0, the CW-WLS smoothing 
method is briefly reviewed and the criteria which serve to the definition of 
data adapted radial parameters of the smoothing windows are discussed. The 
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description of the median smoothing is done in Section An illustration on 
experimental data is reported in Section f|. 

2 CW-WLS Statistical Noise Smoothing with 
Optimal Radial Parameters 

2.1 Local smoothing windows 

Let H = (hij) denote the raw LF 2D-ACAR histogram of interest, corrected 
for finite detector aperture and local variations of the detector sensibility || . 
The statistically relevant information within H sharply decreases when one 
goes from the histogram center towards its borders. This feature is preserved 
in the CF histogram H = (hki), obtained from H under proper discretization 
in the Tp x p y p z frame [Q. 

As a consequence, the histogram bins can be divided into two classes: 
central and border bins. The separation line between the two bin classes is 
somewhat arbitrary. For instance, if we decide to restrict the analysis to the 
projection of the first Brillouin zone onto the histogram plane (1BZ), then the 
central bin area has to include both 1BZ and a surrounding neighborhood of it, 
half 1BZ wide say, to get detailed information on the momentum distribution 
around both sides of the 1BZ boundaries. If, however, the decision is taken to 
include the maximum possible histogram area into analysis, an appropriate 
definition of the central bin manifold includes the largest possible histogram 
area related by a similarity transform, of factor (2A + 1) x (2A + 1) to the 
Dq area of 1BZ. (Typically, A equals two or three, depending on the setup.) 

In ||, a least squares fit of noise- free data to the noisy H data was 
proposed to be performed, at the fractionary coordinates (£, rf) of interest 
inside each (K, L)-th central bin (where the quantities £ and n may vary from 
bin to bin), by means of local approximating surfaces, Ckl, characterized by 
the following two features: 

1. Each Ckl is to consist of an integer number of bins. (This requirement 
follows from the data discretization into square bins, that is, into finite 
regions inside which the structural details of the momentum distribu- 
tion have been averaged out — and hence lost — within the process 
of data acquisition. The bin size is taken henceforth for the unit 
length of the distances in the histogram plane.) 
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2. In the limit of an in-plane continuous point distribution, the approxi- 
mating surface Ckl is to be circular around the point of coordinates 
(K,L). (Such a choice shows the largest linear dimensions under the 
smallest area and it fits isotropically the various possible neighbour- 
hoods. It is thus expected to secure, among the possible 2D shapes, the 
least L 2 norm departure from arbitrary input surfaces.) 

As a consequence, to perform statistical noise smoothing, we draw around 
each (K,L)-th central bin a smoothing window Ckl of quasi-circular shape 
which includes inside it all the bins the centres («, A) of which satisfy the 
inequality 

(k — K) 2 + (A - L) 2 < (2r + l) 2 /4, (1) 
where the quantity r denotes the window radial parameter. 



2.2 CW-WLS smoothing formula 

To accommodate both the data discretization into bins and the possibility 
to predict noise-free values inside the bins, the approximating space of noise- 
free data is spanned by a basis set of polynomials of continuous variables, 
P m (x,y), orthogonal over Ckl- In what follows, it is assumed that the basis 
polynomials are of at most 3-rd degree in each of the variables x and y. 

Over a local window Ckl, the standard deviations aK+k,L+i associated to 
the elements h K +k,L+i , {K + k,L + I) G Ckl of H show little variation 
from (Tkl-i the standard deviation of Ukl at the center of the window ||, [13 



We may therefore assume constant bin weights, (TK+k,L+i = <?kl, within each 
smoothing window (|l|), a hypothesis which results in a constant weight WLS 
smoothing formula of radial parameter r (CW-WLS(r)). 

At the fractionary coordinates (£, if) inside the (K, L)-th bin, 

-0.5 < f < 0.5, -0.5 < 7] < 0.5, (2) 

the CW-WLS (r) yields a smoothed value 

Sk+^l+t, = h K +k,L+i G(£,ri;k,l), (3) 

K+k,L+l£C KL 

where the Green matrix of the smoothing formula is given by 

M 

G&mkJ) = J2 Jom(kJ) P m (^v)- (4) 

m=0 
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Here, u m denotes the norm of the m-th polynomial P m over its definition area 
Ckl, while J mm >(k,l), the overlap integral of the basis polynomials P m (x,y) 
and P m '(x,y) over the area of the (k, Z)-th relative bin inside the current 
smoothing window. For a set of basis polynomials of at most 3-rd degree, 
the upper summation index in Eq. (jij) is M — 10. 



2.3 Optimal values of the window radial parameters 

To achieve the best data fit over an admissible class of functions, a consis- 
tent least squares fit procedure involves two kinds of free parameters: (i) the 
fit parameters which secure minimum L 2 norm departure of the data from a 
particular function belonging to the admissible class, and (ii) the free param- 
eter which selects the best function within the admissible class. The difficult 
point is to fix the last parameter of the procedure. 

The CW-WLS procedure for statistical noise smoothing described in and 
summarized above solves the problem (i) (i.e., it defines the best parameters 
of the constant polynomial degree smoothing filter ([|) under a given radial 
parameter, r, Eq. ([J)). The solution of problem (ii), which consists in the 
definition of the best r value, that is of the most suitable smoothing window 
Ckl under an existing neighborhood of the reference bin (K,L), will be 
discussed below. 

Within the usual polynomial fit, which adjusts a polynomial function of 
variable degree to an input set consisting of a fixed number of elements, to 
get the best polynomial degree, the method starts with a low polynomial 
degree which is then gradually increased by unity until a stopping criterion 
is satisfied. 

Of interest for the present investigation is the criterion proposed by Ham- 



ming |T|j. Considering one-dimensional data, Hamming recommends the 
definition of the best polynomial degree from the study of the distribution of 
the signs of the residuals of the smoothed data. When the polynomial degree 
is increased, a threshold value is reached at which the signs of the residu- 
als over the set under study show stochastic distribution (i.e., the residuals 
mainly consist of statistical noise), while at higher polynomial degrees, the 
prevalence of a given sign over the set of residuals is obtained (that is, a sig- 
nificant part of the useful signal is chopped by the smoothing polynomial). 
The threshold defines the degree of the best fit polynomial. 

While the problem considered by Hamming is different from the present 
one, its principle is very well suited to the definition of the radial parameter, 
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r, Eq. (|T|). Thus, instead of varying the degree of the smoothing surface, 
we adjust its extension by variable r values in a range {r min ,r max }. We 
choose r m i n = 2, which is the lowest available radial parameter value for the 
window (|1|), while r max = 13, a value which corresponds to the chopping of 
the most part of the useful signal by the smoothing procedure. A threshold 
value results, r opt e {r m i n ,r max }, with similar properties to those of the best 
polynomial degree within Hamming's procedure. 

The residual of the smoothed value at the fractionary (£, rj) coordi- 
nates inside the (K + k,L + l)-th bin belonging to the window Ckl, Eq. ([[J) 
is, by definition, 

$K+k+Z,L+l+v = S K+k+^L+l+ri ~ hx+k,L+l- (5) 

Given a Ckl window, signs of quantities (0) are computed over two- 
dimensional sub manifolds of it which obey to two requirements: the sub- 
manifolds have to be large enough such as to result in statistically relevant 
information on the distribution of the signs of the residuals around the ref- 
erence bin, while low enough to secure a sufficiently fast algorithm. 

The least possible manifold of interest is Vi(K,L), the neighborhood of 
the (K, L)-th bin which consists of the reference bin and its eight nearest 
neighbors. If all the nine signs entering Vi(K,L) are identical with each 
other, then the distribution is assumed to be non-stochastic. 

If this rule fails, then we consider the larger neighborhood V^iK^L) of 
the (K, L)-th bin, defined by the 21 bins of the Ckl window (|T|) at the radial 
parameter value r = r min = 2. If two-thirds at least of the signs of the 
residuals inside V2(K,L) are identical, then the statistics is assumed to be 
non-stochastic. 

These two sign count criteria following from Hamming's procedure, how- 
ever, do not exhaust the manifold of non-stochastic two-dimensional distribu- 
tions. There is still the possibility of ordered sign distributions characterized 
by nearly equal occurrence of positive and negative signs. To pick them out, 
the analysis of directional sign distribution of the residuals inside V2(K, L) is 
performed along one or several of the following four directions: the p x axis, 
the p y axis, the first bisectrix, and the second bisectrix. In each case, we 
count the ordered triplets, i.e., the sets of identical three signs at neighbor- 
ing bins along rows, columns, or diagonals respectively. If the count yields 
a majority of ordered triplets with respect to the total number of possible 
ordered triplets in the considered direction, then we decide that the sign 
distribution is non-stochastic. 
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An accidental fulfillment of one of the three abovementioned criteria is 
possible, with the consequence that a spurious r opt cutoff is obtained. To 
avoid such a premature end of the analysis, the requirement of preservation 
of the non-stochasticity of the sign distribution within Vi(K,L) or Vi{K,U) 
under further increase of the radial parameter r is imposed. 

Within an ideal, perfect single crystal, the occurrence of jumps in the op- 
timal radial parameters closely follows the characteristic Fermi surface jumps 
in the measured momentum distribution. Unfortunately, the real single crys- 
tals show a great many number of imperfections of various kinds, which sig- 
nificantly alter the obtained optimal window radial parameters. The crystal 
imperfections blur the Fermi surface jump pattern and make it indistinguish- 
able from a stochastic momentum distribution. Thus, in real crystals, the 
abovementioned definition of the optimal window parameters is to be supple- 
mented with a procedure able to restore the jumps at contiguous bins which 
are characteristic to the Fermi surface. 

Such a procedure is discussed in the next section. 



3 Elimination of Spurious Impurity Effects 
by Median Smoothing 



The median smoothing is known to preserve the jumps, if any, within the 
distribution under consideration. Consistent median smoothing of physical 
quantities associated to an n- axis-projected histogram, H, is obtained pro- 
vided the following two distinct problems are correctly solved. The first con- 
cerns the parameters of the two-dimensional smoothing window: the exten- 
sion of the neighborhood of the (K, L)-th bin over which the median smooth- 
ing is to be done and the weights assigned to the bins entering this neighbor- 
hood. The second concerns the symmetry induced relationships among the 
bins of the histogram. Each of these problems is considered in detail below. 

The extension of the neighborhood entering the median smoothing of the 
(K, L)-th bin can be chosen to be either a Vi(K, L) or a V2(K, L) neighbor- 
hood (defined in Sec. p73| above). 

Within a V\ (K, L) neighborhood, it is natural to assume equal weights to 
all the bins entering it: 
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Within an V^if, L) neighborhood, we can either assume equal bin weights 
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or we can assume different weights. Here we consider the hat shape 
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In the case of data characterized by a high level of crystal imperfec- 
tions, the neighborhood Vi(K,L) with equal weights (f|) is too small to re- 
sult in effective cut of the fluctuations residing in the imperfections. The 
neighborhood V 2 (K,L) with equal weights (|7|) is also inappropriate since it 
overemphasizes the far staying bins, while the weight of the reference bin is 
negligibly small. The most adequate seems to be the choice V2(K,L) with 
unequal weights (§) which includes sufficient enough neighbors while secur- 
ing a fifty-fifty ratio between the weights assigned to the reference bin and 
its nearest neighbours on one side and the farther staying bins on the other 
side. 

Since the information stored in different bins comes from a crystal, there 
are strong symmetry induced correlations of the Fermi-surface-related infor- 
mation stored in a 2D-ACAR histogram. As a consequence, our smoothing 
method has to be able to emphasize it. To implement this aspect into the 
median smoothing, we make reference to the Neumann principle (see, e.g., 
|n>|), according to which, the symmetry originating in positron annihilations 
with electrons within bands crossing the Fermi level, is a subduction to the 
histogram plane of the symmetry point group of the reciprocal lattice of the 
crystal. 

For an n- axis-projected 2D-ACAR histogram of interest here, the distinct 
symmetry elements characterizing this electron fraction are: the inversion 
symmetry center V (placed at the zero-momentum projection of the distri- 
bution), a rotation axis along n, and in-plane symmetry lines. 
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The occurrence of the symmetry classifies the various bins entering an 
n-axis-projected CF histogram H into stars of symmetry. In the case of 
an YBaiCu^O-j^t, single crystal for instance, the stars of the bins entering 
both the c-axis-projected and a-axis-projected histograms contain one, two, 
or four elements, according to the fact that the bin of interest contains the 
symmetry inversion center T inside it, or it lies along a symmetry axis, or it 
is a general bin. 

In conclusion, there are three essential steps which secure the correct 
approach to the median smoothing of an n-axis-projected 2D-ACAR histogram: 

• Derivation of the CF histogram H from the LF histogram H; 

• Derivation of the window optimal radial parameters; 

• Median smoothing over manifolds of data consisting of symmetry stars 
of bins, within neighborhoods V2, with bin weights 

Such a procedure finally yields window optimal radial parameters char- 
acteristic to a structure showing the n-axis-projection of the Fermi surface of 
the sample. As a consequence, if the accumulated statistics contains charac- 
teristic Fermi surface jumps indeed, then the procedure should be able to put 
them into evidence as characteristic jumps in the final window optimal radial 
parameters. That is, while a Fermi-surfaceless structure is characterized by 
gradual variations of the window radial parameters (i.e., Sr jumps between 
neighboring bins in the range {-1, 0, 1}), the occurrence of lines of Fermi 
surface within the 2D momentum projection should be evidenced by jumps 
\5r\ > 2 at contiguous bin positions. 

The next section illustrates this feature on a set of data open to contra- 
dictory interpretation within the usual off-line processing methods. 



4 Results and Their Discussion 

The use of the present method is illustrated on a c-axis-projected 2D-ACAR 
histogram labelled "580", measured on an Y Ba^Cu^O-j-^ single crystal, at 
room temperature, at the University of Geneva, with an upgraded version of 
the setup described in |16j and previously analyzed with different methods 
of resolving the anisotropic component of the 2D-ACAR spectrum || [12|, 17, 
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The parameters which define the crystal frame Tp x p y p z within the labo- 
ratory frame Op^p^pf are 

{% = + fo, 7y = A + Vo, 00, , V'o}- (9) 

Here (7s, Ty) denote the in-plane T coordinates, while <fio, 9 , ip are the 
Euler angles which define the rotations from the Op^PyP® frame axes to the 
Tp x p y p z frame axes. 

The rough approximations (ko,Xo), to the coordinates of T have been 
found to be (—1.0,5.0), while the fractionary coordinates were found to be 
£ = -0.1146 x 10" 1 and r] = 0.2413 (in a D units). The value of the Euler 
angle #o was found to be negligibly small (9 = —0.316 x 10~ 3 radians), 
such that the experimental histogram can be accepted for projected 
histogram indeed. Then the only relevant Euler angle is the sum of the angles 
(f) Q and ^o, «o — —0.12297 x 10~ 2 radians. 

The final result of the procedure described in the previous sections is 
shown in figure |J Both the electron ridge crossing the first Brillouin zone 
from r to X and the pillbox around the S corners are resolved. It is to 
be stressed that the ridge-jumps have been of high amplitude (5r inbetween 
three to six units large), while the pillbox-jumps have always been charac- 
terized by the minimal 5r = 2 value, telling us that the collected ridge jump 
signal is strong, while the pillbox jump signal is weak. 

Whereas the occurrence of the Fermi surface ridge in the analyzed " 580" 
histogram was confirmed by different signal processing techniques in all the 
five references mentioned above, || [T2], |T7], ITS], IS], the only references which 
analyze the pillbox are |19[] and fl~2|| . 



The methods of analysis used in JT9|) did not confirm the occurrence 



of the pillbox in YBa2CusOjs at room temperature. The measurements 
of the positron mean free path in this compound at different temperatures 
suggested that, at room temperature, there is a shallow positron trapping 
mechanism which is responsible of the absence of the pillbox. 

Our previous analysis of the same spectrum [|I2|] has found that the re- 
solved pillbox exists, but it is small (the depth of the Lock-Crisp- West (LCW) 
folded [p0| signal at the S point was found to be only four times larger than 



the experimental errors). 

The present investigation, which is not symmetric component subtrac- 



tion sensitive, adds considerable weight to the conclusions of reference |12 



without making use of the LCW folding to achieve signal enhancement. It lets 
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Figure 1: Occurrence of jumps of the window optimal radial parameters of 
the CW-WLS method for statistical noise smoothing resolves two Fermi surface 
sheets: the ridge, going from r to X across the whole first Brillouin zone of 
the sample and the pillboxes, around the corners S of the first Brillouin zone. 



us infer that, while weakening considerably the pillbox signature at room 
temperature indeed, the positron shallow trapping does not rule it out alto- 
gether. It is to be noted, however, that, although describing a closed contour 
jump in the measured positron-electron momentum distribution, the shape 
of the resolved pillbox around the S point is irregular, a feature which tells 
us about is small amplitude, which can be easily distorted by the occurrence 
of non-homogeneously distributed crystal defects. 
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